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We discuss the QCD phase diagram in the strong coupling limit of lattice QCD by using a new 
type of mean field coming from the next-to-leading order of the large dimensional expansion. The 
QCD phase diagram in the strong coupling limit recently obtained by using the monomer-dimer- 
polymer (MDP) algorithm has some differences in the phase boundary shape from that in the mean 
field results. As one of the origin to explain the difference, we consider another type of auxiliary 
field, which corresponds to the point-splitting mesonic composite. Fermion determinant with this 
mean field under the anti-periodic boundary condition gives rise to a term which interpolates the 
effective potentials in the previously proposed zero and finite temperature mean field treatments. 
While the shift of the transition temperature at zero chemical potential is in the desirable direction 
and the phase boundary shape is improved, we find that the effects are too large to be compatible 
with the MDP simulation results. 
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1. Introduction 



The strong coupling lattice QCD (SC-LQCD) has been successful from the beginning of the 
lattice QCD formulation. In pure Yang-Mills theory, the strong coupling limit (SCL) expression 
of the string tension was proven to follow the area law, and the Monte-Carlo simulation results 
on the string tension were qualitatively explained in the strong coupling region by including var- 
ious plaquette configurations [jj]]. Recently SC-LQCD for pure Yang-Mills theory is extended to 
finite temperatures [ph. SC-LQCD is also powerful in describing the QCD phase diagram. The 
effective potential including fermion contributions were obtained as a function of the chiral con- 
densate in SCL, and the spontaneous chiral symmetry breaking in vacuum and its restoration at 
finite temperature and density have been discussed based on the derived effective potential [^]- 
|6p. Recently, SC-LQCD framework with fermions is extended to include the next-to-leading order 
(NLO, ff(\/g 2 )) [^,|§] and the next-to-next-to-leading order (NNLO, @{\jg A )) @ contributions 
of the strong coupling expansion. The mean field approach in SC-LQCD overestimates the MC re- 
sults of the transition temperature (T c ) by around 10 % in SCL, while we overestimate T c by about 
50-60 % at /$ = 2N c /g 2 ~ 4. When we include the Polyakov loop effects in SC-LQCD, it is possi- 
ble to roughly reproduce T c or the critical coupling (/3 C ) at a given temporal lattice size (N z = 1/7) 
in the coupling region j8 = 2N c /g 2 < 4 [|l(i|]. 

Having these successes in mind, we expect that the phase boundary predicted in the mean 
field treatment of SC-LQCD is promising in the strong coupling region. Recently de Forcrand and 
Fromm [11 ] obtained the phase diagram in SCL-LQCD in a Monte-Carlo simulation based on the 
monomer-dimer-polymer (MDP) algorithm [12]. In the MDP simulation, we first integrate out the 
link variables in SCL (so-called zero T treatment), and the integral over quarks is simulated by the 
sum over loop configurations. Since the gauge integral is carried out analytically, the sign problem 
is weakened. The shape of the phase boundary is somewhat different from that in the mean field 
predictions; T c for a given pL is much lower in the MDP results. It is important to understand the 
origin of this deviation, since the MDP simulation is only applicable to SCL ( 1 /g 2 = 0) and the 
mean field treatment of SC-LQCD is one of the few approaches in which we can discuss cold dense 
matter directly based on non-perturbative QCD. 

The deviation suggests that some of the approximations adopted in the mean field treatment 
of SC-LQCD would not be good enough at finite jit. There are two types of approximations in 
the mean field approach. One of them is the assumption that the auxiliary chiral field takes a 
constant value, and the other is the truncation in the large dimensional (l/d) expansion [13]. In 
this proceedings, we discuss the possibility to introduce another type of mean field from the NLO 
term of the l/d expansion than the chiral condensate in SCL. Specifically, we examine the role of 
the mean field of the type V± VtX = f]v,xXxXx+v in the zero T treatment. 



2. Strong coupling limit of lattice QCD with another mean field 

In this section, we first review the framework of the mean field treatment of SCL-LQCD 
briefly, and compare the obtained phase diagram with that in the MDP simulation. Next we in- 
troduce another kind of mean field. Throughout this proceedings, we set the lattice spacing as 
a = 1. 
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We consider the lattice QCD action with one species of unrooted staggered fermion for color 

SU(JV C = 3), 



iT LQ cD=/ @[X,X,U V ] e- SF - Sa , 



(2.1) 
(2.2) 



x v=0 



where x(x)> m 0> U vx , r\ VJC = exp(5 v o/x)(— \) Xo+ - +x i-^ denote the quark (antiquark) field, the bare 
quark mass, link variable, and the staggered phase factor, respectively. The pure Yang-Mills action 
So is proportional to 1 /g 2 , and disappears in SCL, g — > oo. 
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Figure 1: Phase boundaries in the MDP simulation (dashed lines) Jl 1|], the standard mean field treatment 
(dotted lines) [^], and the present treatment (solid lines). Left (right) panel shows the phase boundary in the 
lattice unit (normalized by T c (^i = 0) and fj, c (T = 0)). 

In the finite T treatment of SCL-LQCD, referred to as the standard mean field MF (Std.) in 
the later discussion, we first integrate out spatial link variables, and obtain the effective action of 
quarks and the temporal links [Eh, 



n(std.; 

>eff 



1 



1 



^1 e^XxUo,xX x+ 6-^X x+ dUl x Xx ~ jrr £M,M* + s +m £M v + 0(1 /Vd) , (2.3) 



where d = 3 is the spatial dimension, and M x = XxXx represents a mesonic composite. We next in- 
troduce an auxiliary field for the chiral condensate to make the effective action bilinear in fermions, 
and integrate fermions and temporal link variables, for which the anti-periodic and periodic bound- 
ary conditions are imposed, respectively. The effective potential as a function of the chiral auxiliary 
field a is obtained in the mean field approximation as, 



9. 



eff 



4N r 



riog 



&wh{(N c + l)E q /T) 
sinh(E q /T) 



+ 2cosh(N c lx/T) 



(2.4) 



where E q = arcsinh(<ia/2A^ c + mo) is the one-dimensional quark excitation energy. The phase 
diagram in MF (Std.) is shown by the dotted curve in Fig. |TJ Compared with the MDP results, 
we overestimate the transition temperature at /I = by 10-20 %, while we roughly reproduce 
the critical chemical potential at T = 0. In between, the transition temperature stays high in the 
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chemical potential region /I < 0.5 in MF (Std.), and the phase boundary shape is significantly 
different from that in MDP. 

There are two types of approximations adopted in MF (Std.). One of them is the assumption 
that the introduced auxiliary field a is constant, and the other one is the truncation in the l/d 
expansion. In the l/d expansion [13], the M x M x+ y term containing four quarks and sum over di- 
mensions is assumed to be finite at large d. Then the quark field scales as d~ l l A , and the NLO 
terms containing six quarks are proportional to 1 j\fd. Baryonic composite action belongs to NLO, 
and it is natural to expect that the baryonic action would affect the phase boundary at finite \i. 
Baryonic composite effects have been discussed at finite T in Refs. [14], but the adopted bosoniza- 
tion method [ p3| ] is not fully compatible with the chiral symmetry. Further studies are necessary to 
discuss the chiral phase transition seriously. 

By comparison with MF (Std.), both spatial and temporal link variables are integrated out first 
in the zero T treatment. The effective action including the NLO terms in the l/d expansion is given 
as. 



>eff : 



— ^M X M X+V + m £M V + - £ [nl x B x B, 



r\-% + M+e{lHd + l)), (2.5) 



x+v ' 



where B x = £ a bcXxXx Xx/^ * s tne baryonic composite. The third term in Eq. (2^1 is the NLO 
term, @{l/ \fd), in the l/d expansion; it contains six fermions (oc <i~ 3 / 2 ) and sum over space-time 
dimensions. This effective action with higher order terms in the large dimensional expansion is 
used in the MDP simulation [12]. In Ref. rf|], the phase transition at finite ji. and zero T was 
investigated by introducing the auxiliary baryon field b ~ B in the mean field approximation for 
the chiral condensate, but the phase transition at finite T is not well described. This shortcoming 
could come from the lack of the anti-periodic boundary condition, which is decisive in the chiral 
transition at finite T. 

In order to impose the anti-periodic boundary condition for quarks in the effective action 
Eq. (2.5), we here consider another type of auxiliary field, which corresponds to the point-splitting 
mesonic composite, 



V. 



+ V.X 



^Iv.xXxX. 



,X+V ) 



T]v,xXx+vXx 



(2-6) 



We note that the baryonic term in Eq. ( |2.5| ) can be rewritten by using the anti-commuting property 
of the Grassmann variables as 



Vv^BxBjc 



(V +v , x ?/6 



v.x^x+v 



By 



(V- V ,x?/(> 



By applying the extended Hubbard-Stratonovich transformation 



exp(=faV 3 ) ~exp 
(i)y2 



exp(ai/A 



^exp 



-a(V/ (2 V (2) - 



.yV (1) ±V^ (1 V) 



we obtain the effective action for quarks and auxiliary fields, 



Z± v =2a ( vdv =F Y±l =F Y±lv±l 



(2.7) 

(2-8) 
(2-9) 

(2-10) 
(2.11) 
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Under the assumption that auxiliary fields take constant values, we can evaluate the Matsubara 

Ik) 

product with the anti-periodic boundary condition for quarks. Equilibrium condition for y±v an d 

- (k) 

Y±v is used to reduce the number of independent variables. For example, equilibrium values of 



Y+l and are related as, \}f+v 
rotational and reflection symmetry, <p + j 
be as follows. 



q> +v and y/ 



(i) 



(2) 

-Y+v 



We also utilize the 



= (p s (j = 1,2,3). The effective potential is found to 



eff ; 



d Jrl c 2 + 2a{(pl + (pl)+4ad(p i s +V q (a = 1/48) , 



^ + log (l +e-^-^ T ) +log (l + e -(^W/r 



Z± =6aq>± , 
Ek =arcsinh(£jt/Z z ) , Ek 



Z x = ^Z + Z_ , fl = H + log(Z+/Z_ 



2 + Z2sin 2 k 



in 



m„ 



d + l 

~2N r 



-G + rriQ 



N c \ogZ x , 



(2.12) 
(2.13) 
(2.14) 



6a <p; . (2.15) 



This effective potential has interesting features; quarks couple with auxiliary fields via the con- 
stituent quark mass m q and the wave function renormalization factor Z±, s , and it contains the mo- 
mentum integral with a lattice type replacement, k — > sin/:. 

It should be noted that the composite V± v is not gauge invariant, and the mean field introduced 
here should be regarded as the one in a fixed gauge. We expect that this feature may not be serious, 



since we are discussing the dynamics in the link integrated effective action Eq. (2.5) 



3. Effective potential surface and phase boundary 



We shall now examine the effective potential Eq. ( [2.1 5[ ). For simplicity, we ignore the effects 
of (p s which connects spatially separated quarks, and we discuss the results only in the chiral limit, 
niQ = 0. We first consider the j£" e ff in vacuum, (T,/j.) = (0,0). At pL = 0, the temporal forward 
and backward auxiliary fields take the same equilibrium value, q> + = <p_. When (p± are zero, the 
vacuum effective potential becomes the logarithmic type, J^^f — > b a o 2 /2 — N c \ogm q , which is the 
leading order effective potential of the \/d expansion in the zero T treatment [j||. For finite <p±, 
the effective potential becomes the arcsinh type, which is typical in the finite T treatment. Thus the 
present effective potential may be regarded as an interpolating one of the zero and finite T effective 
potentials. 

In Fig. ^[ we show the effective potential as a function of a and cp = ^/cp + (p- at several 
(r,/x). In vacuum, equilibrium is realized at finite a and q>± = 0, while <p± (a) grows (decreases) 
as T increases. The phase transition at finite T and zero /I is the second order. The transition 
temperature at \i = (T c = 0.92) is smaller than those in MDP (T c ~ 1.4) and MF (Std.) (T c = 
5/3). The decrease of T c may be understood as the contribution from single quarks. In MF (Std), 
coherently moving three quarks contribute to the effective potential as a baryon as seen in the 



Boltzmann factor of exp(— N c /J./T) in Eq. (2A). By comparison, the mean field (p± allows a single 



quark excitation as found in the Boltzmann factor exp(— ju/r) in Eq. ( 2.13 ). 

At finite /I and T = 0, the energy surface is separated by the ridge at E q = p., and the vacuum 
configuration (finite a and zero q>) jumps to the high density configuration (zero a and finite cp). 
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r eff 




Figure 2: Effective potential surface in the (a, (p = ^/<p + <p_) plane. The solid line connects the points which 
give the lowest J£" e ff for a given <p value, and the filled circles show the equilibrium. 



The transition chemical potential at T = (;ii c = 1.08) is larger than those in MDP (jj, c ~ 0.59) and 
MF (Std.) (ix c = 0.55). 

In Fig. [j], we compare the phase boundary in the present treatment including (p± with those 
in the MDP simulation and MF (Std.). The new mean fields q>± shift T c and pL c from those in 
MF (Std.) in right directions, and when we normalize T and /I by r c (ju = 0) and jU c (r = 0), 
the phase boundary shape is significantly improved and becomes is similar to that in the finite T 



treatment with baryonic composite effects Q14|]. However, their effects are too much in the shifts of 
T c and jX c . 



4. Summary 



In this proceedings, we have discussed how we can understand the phase boundary in the 
Monomer-Dimer-Polymer (MDP) simulation in the strong coupling limit (SCL) of lattice 
QCD for color SU(3) with unrooted staggered fermion in an analytical method based on the mean 
field treatment. Since the MDP simulation is available at present only in SCL, it is important to 
understand it in a method which is applicable to finite coupling cases. Here we have examined a 
new type of mean field <p±, which connects the fermion in a different temporal variable in the zero 
T treatment, where both spatial and temporal link integrals are carried out first. This enables us to 
impose the anti-periodic boundary condition of quarks, and may be important to describe the phase 
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transition at finite T. Actually we can describe the finite T phase transition in the zero T treatment. 
The new mean field shifts the transition temperature at /i = and the transition chemical potential 
at T = in the right directions, but we find that the effects are too much to be compatible with 
the MDP results. It is still a challenge to understand the phase diagram in the MDP simulation. It 
would be necessary to examine the combination with the baryonic effects proposed so far []|, [l4[] 
and the fluctuation of the chiral condensate. 
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